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We study representations of the Schrodinger algebra in terms of operators in non- 
relativistic conformal field theories. We prove a correspondence between primary 
operators and eigenstates of few-body systems in a harmonic potential. Using the 
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correspondence we compute analytically the energy of fermions at unitarity in a har- 
monic potential near two and four spatial dimensions. We also compute the energy 



of anyons in a harmonic potential near the bosonic and fermionic limits. 
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I. INTRODUCTION 

Conformal field theories (CFTs) form a special class of relativistic quantum field theories, 
where the Poincare symmetry group is enlarged to the group of conformal transformations. 
One element of the conformal algebra is dilatation: CFTs are always scale invariant. The 



> . 

conformal algebra and its representation have been studied extensively [1]. 
! In this paper we study nonrelativistic counterparts of relativistic conformal field theo- 

ries. There are several examples of such theories beside the trivial noninteracting theories. 
Nonrelativistic particles interacting through a 1/r 2 potential is one such example. The phys- 



ically most important example in three spatial dimensions is the theory describing spin-1/2 
fermions with point-like interaction fine-tuned to infinite scattering length (fermions at uni- 



tarity) [2]. Such fermionic systems have been created and studied experimentally. Theories 
describing anyons present another class of nonrelativistic CFTs, but in two spatial dimen- 
sions. 

The nonrelativistic analog of the conformal algebra is the so-called Schrodinger alge- 
bra [3, 4]. While the Schrodinger algebra has been studied before [5-7], we are mostly 
interested in the representations of the Schrodinger algebra in terms of operators. We show 
that the concept of primary operators can be directly transferred to nonrelativistic theories. 
In addition, we show that there is an operator-state correspondence: a primary operator 
(with some exceptions) corresponds to an eigenstate of a few-particle system in a harmonic 
potential. The scaling dimension of the operator coincides with the energy of the corre- 
sponding eigenstate, divided by the oscillator frequency. 

The operator-state correspondence allows us to translate the problem of finding the en- 
ergy eigenvalues of a few-fermion state at unitarity, or a few-anyon state, in a harmonic 
potential to another problem of finding the anomalous dimensions of composite operators 
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in the nonrelativistic conformal field theory in free space. The latter problem is amenable 
to standard diagrammatic techniques for fermions at unitarity near two or four spatial di- 
mensions, or for anyons near the bosonic and fermionic limits. We present a few examples 
of such calculations in this paper. In particular, we compute the ground state energy of 
up to six fermions at unitarity in a harmonic potential near two and four dimensions, and 
interpolate the results to find the energy in three dimensions. We also compute the ground 
state energy of up to four anyons in a harmonic potential. 



II. SCHRODINGER ALGEBRA 



A. Derivation of the algebra 

We briefly review the Schrodinger algebra [3, 4]. For definiteness, consider a nonrelativis- 
tic theory described by a second-quantized field ip a (x) (where a is the spin index) which 
satisfies the commutation or anticommutation relation 

bl> a (x), il>l(y)}± = 6(x - y)5 af3 . (1) 

Throughout this paper, we use nonrelativistic natural units h — m — 1 where m is a particle 
mass. We consider a general spatial dimension d. Define the number density and momentum 
density, 

n{x) = ^(x)i>(x), ji(x) = - l -(^(x)d^(x) - d^{x)^{x)) (2) 
(summation over spin indices is implied). Their commutators are 

[n(x), n(y)} = 0, [n(x), ji(y)} = -in(y)di5(x - y), (3a) 
[ji(x), jj(y)} = -i (jj{x)di + ji(y)dj) S(x-y). (3b) 

The Schrodinger algebra is formed by the following operators: 

N = Jdxn(x), Pi = Jdxji(x), = J dx (xijj(x) - Xjji(x)), (4) 

Ki = Jdxxin(x), C— Jdx^-n(x), D = Jdxxiji(x), (5) 

and the Hamiltonian H . The operators in Eq. (4) have simple physical interpretation: iV 
is the particle number, Pi is the momentum, and is the orbital angular momentum. In 
a scale-invariant theory like unitary fermions, these operators form a closed algebra. All 
commutators except those that involve H can be found from Eqs. (3). First iV commutes 
with all other operators: 

[N, any] = 0. (6) 
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The commutator of the angular momentum M^ with an operator is determined by the 
transformation properties of the latter under rotations, 

[M tj , N] = [Mij, C] = [Mij, D] = 0, (7) 
[Mij, P k ] = i(5 ik P 3 - 5 jk P t ), [M^, K k ] = i(5 lk K j - 5 3k Ki), (8) 
[Ma, M kl ] = i(5 lk M jk - 8 jk M a + 5 u M kj - 5 3l M ki ). (9) 

The remaining commutators are 

[Pi, Pj\ = [K h Kj\ = [K h C\ = 0, [K h Pj\ = td t ,N, (10) 
[D,Pi\=iPi, [D,Ki\ = -iKi, [D, C\ = —2iC. (11) 

Now let us find the commutators of the Hamiltonian H with other operators. Conservation 
of particle number, momentum, and angular momentum implies that 

[H, N] = [H, Pi] = [H, M^ = 0. (12) 

On the other hand, the continuity equation implies that 

[H, n] = -id t n = idji, (13) 

from which it follows that 

[H, Ki\ = -iP h [H, C\ = -iD. (14) 

The computation of the last commutator [H, D] requires the condition of scale invariance. 
For definiteness, let us consider fermions at unitarity, described by the Hamiltonian 

H = Jdx l -d^ a d^ a + l - Jdxdyi/>i(x)if>l(y)V(\x-y\)Mv)Mx), (15) 

where V(\x — y\) is a short-range potential with infinite scattering length. We note that D 
is an operator of dilatation, 

e- lXD ^j(x)e lXD = e dX/2 ^j(e x x), (16) 

from which one finds 

e~ lXD He iXD = e 2X H', (17) 
where H' is the same as H but the potential V is replaced with a new potential: 

V(r) -> V'{r) = e- 2X V(e- x r). (18) 

If V corresponds to infinite scattering length, then V also corresponds to infinite scattering 
length. From the point of view of long-distance physics, H' = H. Therefore, we find 

[D, H] = 2iH. (19) 
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It is clear that Eq. (19) is simply the condition of scale invariance, and hence must hold for 
noninteracting anyons in two spatial dimensions. 

A more lengthy proof, which can be given for particles at unitarity, is to use the momen- 
tum conservation equation 

dtji + djUij = 0, (20) 

where Uij is the stress tensor, which can be defined for a generic potential V(\x — y\) (see 
Appendix B). The commutator is then 

[D, H]=iJ dxU u {x). (21) 

At unitarity, one can show that (see Appendix B) 

JdxIUi(x) = 2H, (22) 

and Eq. (19) follows. 

The Schrodinger algebra is summarized in Appendix A. 



B. Local operators and representations of the Schrodinger algebra 

We introduce the notation of local operators 0(t, x) as operators which depend on the 
position in time and space t, x so that 

0(t, X) = e iHt-iPiX i0 ^ e -iHt+iPiXi m (23) 

A local operator O is said to have scaling dimension A if 

[D,O(0)]=tA o O(0), (24) 

and to have particle number N& if 

[N, O(0)} = N o O(0). (25) 

We shall consider only operators with well-defined particle number and scaling dimension. 
Examples of such operators are tp and A^ = A^ = d/2 and — = = 1. In the 
case of fermions at unitarity, a more complicated local operator is 

4>{x) = lim \x - y^-^x^iy). (26) 

V-* x 

The presence of the prefactor \x — y\ d ~ 2 guarantees that the matrix elements of the operator 
4>(x) between two states in the Hilbert space are finite. 1 The scaling dimension of is 



The condition of unitarity requires that the wave function of N spin-up and M spin-down fermions 
^(xi, X2, ■ ■ ■ , xm] 2/i, J/2, • • • , 2/m) behaves like \xi — yj\ 2 ~ d when \xi — yj\ — > for any pair of fermions 
with opposite spins i, j. 
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A^ = 2. This can be found applying elementary dimension counting to Eq. (26): 

= 2A^ + (d - 2)A X = 2^ + (d- 2)(-l) = 2. (27) 

Let us look at the set of all local operators O a (t,x). These operators, put at t — and 
x = 0, form a representation of the Schrodinger algebra: for any operator A in the algebra 

[A, O a (0)\ = A ab O b (0). (28) 

We shall discuss the irreducible representations of the Schrodinger algebra. 

First we notice that if O has dimension A , then [P h O] has dimension A + 1: 

[D, [P U 0\\ = [P U [D, 0\\ + [[D, P^, O] 

= [Pi, iA 0] + [iPi, O] = t(A + l)[Pi, 0\. 

Analogously the dimensions of [K { O], [C, O], and [H, O] are A - 1, A - 2, and A + 2, 
respectively: 

[D, [K i ,0]]=i{A -l)[K i ,0], (30) 
[D, [C,0]]=i(A -2)[C,0], (31) 
[D, [H, 0\\ = i(A + 2)[H, O]. (32) 

Assuming that the dimensions of operators are bounded from below, if one starts with a 
given operator and repeatedly takes its commutator with Ki and C, one lowers the dimension 
of the operator until it cannot be lowered further. The last operator O obtained this way 
has the property 

[Ki, O] = [C, O] = 0. (33) 

Operators that commute with Ki and C will be called primary operators (quasiprimary 
operators in the terminology of Ref. [7]). In general, [Ki, O] — does not imply that 
[C, 0\ — 0, and vice versa. 

Starting with a primary operator O, one can build up a tower of operators by taking 
commutators with Pj and H . In other words, starting with primary operators one can 
build up whole towers of operators by taking their space and time derivatives. For example, 
the operators with dimension Aq + 1 in the tower are [Pi, O] = idiO. At the next level 
(dimension Aq + 2), the following are possible: 

[H, O] = -id t O, [Pi, [Pj, 0]\ = -didjO. (34) 

Commuting the operators in Eq. (34) with Ki and C, we get back the operators in the lower 
rungs of the tower. 

It is easy to see that the operators built from a primary operator by taking space and 
time derivatives form an irreducible representation of the Schrodinger algebra. It is also 
possible to show that the full set of all local operators can be decomposed into irreducible 
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representations, each of which is built upon a single primary operator. The task of finding 
the spectrum of dimensions of all local operators reduces to finding the dimensions of primary 
operators. 

For an operator Q(t, x) with dimension Aq at an arbitrary spacetime point, the following 
commutation relations hold, 

[Pi, O] = idiO, [H, O] = -id t O, (35) 
[D, O] = i(2td t + xA + A )0. (36) 

Moreover, if O is a primary operator then 

[K i ,0] = (-itd i + N x i )0, (37) 

2 

[C, 0\ = -i{t 2 d t + txidi + tA )0 + X -N 0. (38) 

The exponentiated version of Eq. (36) is 

e~ iXD 0(t, x)e iXD = e XA °0(e 2X t, e x x). (39) 

For any set of n operators, one can define an n-point correlation function, 

G n (h, Xl , t 2 , x 2 ,..., t n , x n ) = (0|T(9i(ti, x 1 )0 2 (t 2 , x 2 )... O n (t n , x n )\0), (40) 

where T is time ordering. Clearly for G n to be nonzero it is necessary that No 1 + Nq 2 + 
• • • + No n = 0. If all Oi have definite dimensions then this correlation function has a scaling 
property 

G n (e 2 % e x Xi ) = exp ^-A A o)j G n (t h x t ), (41) 

which follows from Eq. (39) and e iXD \0) = |0). 

The correlation functions of primary operators are further constrained [7] . As an example, 
consider the two-point correlation function of a primary operator O with its Hermitian 
conjugate: 

G(t, x) = (0\TO(t, x)O*(0)\0). (42) 
Using (0\[Ki, TO(x)O*(y)]\0) = and Eqs. (37) one obtains 

(-itdi + N Xi)G(t, x) = 0. (43) 

Combining with the scale invariance, the two-point correlation function is determined up to 
an overall coefficient, 

(44) 



G(t, x) = Cr Ao exp (-iNo^j 
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C. Correspondence to states in a harmonic potential 

We now show that each primary operator corresponds to an energy eigenstate of a system 
in a harmonic potential. We set the oscillator frequency of the harmonic potential w to 1. 
The total Hamiltonian of the system in a harmonic potential is 

H osc = H + C. (45) 

Consider a primary operator O put at t — and x = 0. Let O be constructed from 
annihilation operators, so that acts nontrivially on the vacuum |0). Consider the following 

state 

\m a ) = e - H O j \0). (46) 

If the particle number of is N \, then \$>o) is an iV^t-body state. Let us show that \^o) 
is an eigenstate of the Hamiltonian: H osc = H + C. Indeed 

H osc \*o) = e- H (e H H osc e- H )O*\0). (47) 

We now use the formula 

e H H osc e' H = H osc + [H, H osc ] + ±[H, [H, H osc ]] + ■■■ . (48) 

Using the commutation relations in Appendix A, we find that all terms in the • • • vanish, 
and the right hand side is equal to C — iD. Therefore 

tfosd^o) = e~ H {C - iD)O*\0) = e~ H 0\C - iD)\0) + e - H [C - iD, O*]\0). (49) 

However, both C and D annihilate the vacuum, C|0) = D\0) = 0, and since O is a primary 
operator, [C, C f ] = 0. Thus, using [D , O*] = -[D , O]* = iA a O^, we obtain 

tfoscl^o) = e- H A o O^\0) = A |* ), (50) 

i.e., \^o) is an eigenstate of the system of N t particles in a harmonic potential, with the 
energy eigenvalue A^ (times fiw). 

It is known that the eigenstates of H osc are organized into ladders with spacing between 
steps equal to 2 [8, 9]. The raising and lowering operators within a ladder are [8] 

L + = H -C + iD, (51) 
L_ = H -C -iD. (52) 

Let us show that the state \^o) is annihilated by L_ and hence is the lowest state in its 
ladder. Indeed, using the identity 

e H L_e~ H = -C, (53) 

we find 

L^o) = e- H e H L_e- H O^\0) = - e - H CO*\0) = -e~ H O^C\Q) = 0. (54) 

Clearly, in order to correspond to a nontrivial eigenstate of H osc , must not annihi- 
late the vacuum: C^|0) ^ 0. We shall consider the operators O that are built from the 
fundamental annihilation operators of the field theories. 
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D. Simple examples: one and two-body operators/states 

Let us illustrate this correspondence using one-particle and two-particle operators ip and 
at unitarity. The operator -0 has scaling dimension d/2, which matches the ground state 
energy of one particle in a harmonic potential in spatial dimension d. The operator has 
scaling dimension 2. The ground state of two particles at unitarity in a harmonic potential 
has the wave function 

e -(x 2 +2 / 2 )/2 

(j>(x : y)cx |a ,_ y | d _ 2 , (55) 
and the ground state energy is also 2. 

III. EXAMPLE 1: FERMIONS AT UNITARITY 

In this section, we compute the scaling dimensions of some operators in the theory de- 
scribing spin-1/2 fermions at unitarity. In order to have a small parameter for perturbative 
expansions, we shall work near two and four spatial dimensions, and then, interpolate the 
results to the physical three spatial dimensions. Since the energy eigenvalues of two and 
three fermions in a harmonic potential can be found exactly, we can use these cases to test 
our expansions and interpolation schemes. In the cases of more than three fermions, only 
numerical results exist. Our analytical calculations, as we will see, are consistent with the 
numerical ones. 

There are two field-theoretical representations of fermions at unitarity, one becoming 
weakly coupled as d — > 4 and the other becoming weakly coupled as d — > 2 [10, 11]. We 
shall consider these two cases separately. 

A. Near four spatial dimensions 

1. Fixed point 

In the first representation, the Lagrangian density describing fermions at unitarity is 

C = i^ld^ a - i| W CT | 2 + i<f>*d t <j> - ||V0| 2 + # T Vj0 + gipiM*- (56) 

The canonical dimensions of the fermion field ip a and the boson field are both d/2. There- 
fore the coupling constant g is relevant at weak coupling below d — 4. There are three other 
relevant terms one can add to the Lagrangian density (56): \x a i\)\i\)v and — (g 2 /c ) 0*0. \i a 
is a chemical potential for each spin component of fermions and here we consider the system 
at zero density /i a — 0. Furthermore we assume that the system is fine-tuned so that the 
coefficient in front of 0*0 satisfies 
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-»E(p) 




FIG. 1: One-loop self-energy diagram to renormalize the wave function of 4>. 

[( c o) -1 is zero i n dimensional regularization.] This condition is equivalent to the fine-tuning 
to the infinite scattering length. We denote the propagators of if) and by G(p) and D(p), 
respectively. 

The renormalization of the theory can be performed in the standard way considering 
e = 4 — d to be a small parameter for perturbation. There is a one-loop self-energy diagram 
for which is logarithmically divergent (Fig. 1). Integrating out modes in the momentum 
shell e~ s A < k < A, we obtain 



9 2 ( P 2 \, A 
p - — In ■ 



(58) 



8tt 2 V 4 ; e- s A' 
which corresponds to the wave-function renormalization of 



<f -,s. (59) 



in 2 



The anomalous dimension of is found by the standard formula 

ldln^ g 2 
7 ^ = -2^^ = 16^ (60) 

There is no divergent one-particle irreducible diagram that renormalizes the t/^t^* cou_ 
pling. As a result, the (3 function that governs the running of g, 

| = (61) 

is determined by the dimension of g: 

M=(2- d --„) g =i g -^- 2 . (62) 
There is a fixed point located at 

g 2 = 8vr 2 e. (63) 
At this fixed point, the theory is a nonrelativistic CFT describing fermions at unitarity. 
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FIG. 2: One-loop diagram to renormalize the three- fermion operator tfiipp 

2. Scaling dimensions of operators 
Since the one-fermion operator if> is not renormalized, its scaling dimension is 

= \- (64) 

Using the operator-state correspondence, we find that there is a one-fermion state in a 
harmonic potential with energy (d/2)u, which is obvious. 

The two-fermion operator 0, on the other hand, has a scaling dimension different from 
its canonical dimension d/2. At the fixed point, 

A^ = ^ + 7^ = 2. (65) 

This is, of course, consistent with Eq. (27). Since there is no other contribution to A^, 
Eq. (65) is exact to all order in e. According to the operator-state correspondence, this 
implies the existence of a two-fermion state with zero orbital angular momentum and energy 
2uo. The wave function of this state is given by Eq. (55). 

Let us look at three-fermion operators. The simplest operator is <f)ip-\. This operator 
has zero orbital angular momentum 1 = 0. The diagram that contributes to its anomalous 
dimension to order e is depicted in Fig. 2. It is evaluated as 

/rlk n 2 A 

— G( k)D ( k)G(-k) = -l-^—. (66) 

Therefore the renormalized operator differs from the bare operator by a rescaling factor: 

(M)rcn = Zw,<N>h Where 

g 2 

Z H, = 1 - q^s. (67) 

At the fixed point, —dlaZ^ /ds = 4e/3 is the anomalous dimension of the composite 
operator cf)ip^ (more precisely, the nontrivial part of the anomalous dimension because there 
is a trivial part equal to 7^). We thus find the scaling dimension 



= A^ + A^ T + \e = 4 + jje. (68) 
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According to the operator-state correspondence, (jnjj^ corresponds to a three-fermion state 
with I = and energy equal to 



E (0) - 



4 + - e + 0(e 2 ) 



UJ. 



(69) 



This state continues to the first excited state of three fermions in a harmonic potential at 
d = 3. Even within the leading correction in e, the result extrapolated to e = 1, pa 
4.83 u, is not too far from the true result of 4.66622 c; at d — 3. 

The three-fermion ground state in a harmonic potential at d = 3 has 1 = 1. There are two 
lowest I = 1 operators near four spatial dimensions; 0V^t and (V^)V't- Let us consider the 
renormalization of a general operator acp'Vipi + b (V0). Inserting this operator into Fig. 2, 
we find 



ig 



dk 



{2n) d+1 
a + 5b 



6tt 2 



k + q) + b (k + p)} G(k)D(k + p)G(-A; + q) 
A 
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5a + 76 
P+ ^2- 9 



(70) 



In 



5 A 



In order to have well-defined anomalous dimensions, this should be proportional to ap + bq. 
Thus we have two solutions (a, b) oc (1,1) and (a,b) oc (2,-1), which have anomalous 
dimensions 4e/3 and — e/3, respectively. The first possibility corresponds to the operator 
V(0 , 0f), that is obviously not a primary operator. Its scaling dimension is trivially equal to 
A^ T + 1 (corresponding to an excitation in the center of mass motion). The other operator 
20V , 0t ~ C^ 7 0)V ; t h &s the nontrivial scaling dimension 

1 5 

A20Vi/< t -(v</>)v>t = A<£ + A^ T + 1 — -e = 5 — -e. 

The operator-state correspondence tells us that the three-fermion state with I 
harmonic potential has the energy 



E. 



(i) 



5-^ e + 0(e 2 ) 
6 



UJ. 



(71) 
1 in a 

(72) 



The extrapolation to e = 1 gives E^ pa 4.17 u, which is not too far from the true ground 
state energy 4.27272 ui at d = 3. 

We now turn to the four-fermion state with I = represented by the operator 2 . The 
first nontrivial correction to its scaling dimension is of order e 2 given by the diagram depicted 
in Fig. 3. The two-loop integral can be performed analytically and we find 

27 



Ax> = 4 + 8e z In 



16' 



Thus the ground state of four fermions in a harmonic potential has the energy 



E 



(o) 



27 

4 + 8e 2 ln- + 0(e 3 ) 
16 



UJ. 



(73) 



(74) 
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p + q 



-p-q 



FIG. 3: One-loop diagram to renormalize the four-fermion operator </> 2 . 

The correction, although it is of order e 2 , has the large coefficient. Therefore, in order to 
extrapolate to e = 1, we shall not use Eq. (73) directly but will combine it with a result 
near two spatial dimensions. 

For a general even number of fermions N = 2n, the operator <ft n corresponds to the ground 
state in a harmonic potential with / = 0. Its scaling dimension is given by 

A r =N + N(N-2)e 2 \n^ + 0(e 3 ). (75) 

When N = 2n+l is odd, the operators (p n ip-\ and 20 n V^|— 71 ~ 1 (V0) , ?/>t correspond to states 
with orbital angular momentum I = and 1 = 1, respectively. Their scaling dimensions are 

AN — 7 

A^ T = N + 1 + — — e + 0(e 2 ) (76) 

and 2N-21 

A 2< /,"vV'T-</' n " 1 (v^)V'T = N + 2-\ — — e + 0(e ). (77) 

lo 

According to the operator-state correspondence, the energy of A-fermion state in a harmonic 
potential is simply given by E$ = Ao w. 

The leading-order results [E^ = N uo for even N and E$ = (N + l)u and E$ = 
(N + 2) uj for odd N] can be easily understood by recalling that, in the limit of d — * 4 from 
below, fermion pairs at unitarity form point-like bosons and they do not interact with each 
other or with extra fermions [10, 12]. So the ground state for N = 2n fermions consists of 
n free composite bosons, each of which has the lowest energy 2 uj in a harmonic potential at 
d = 4. When N = 2n + l, the ground state has / = and consists of n composite bosons and 
one extra fermion in the lowest energy states. In order to have an I = 1 state, one of the n+1 
particles has to be excited to the first excited state, which costs additional lu. At d = 4, we 
observe the odd-even staggering in the ground state energy as E$ — (En-i+En + i) /2 = 1u> 
for odd N. 
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FIG. 4: One-loop diagram to renormalize the four-point vertex coupling g 2 . 
B. Near two spatial dimensions 

1. Fixed point 

The other field-theoretical representation of fermions at unitarity is provided by the 
Lagrangian density 

C = i^dtfa - i|Wv| 2 + sVfyfyl^T- (78) 

Again we consider the system at zero density fi a = 0. When e = d — 2 is small, the 
coupling g 2 is renormalized by the logarithmically divergent one-loop diagram (Fig. 4). The 
renormalization group equation for g 2 is 

The fixed point is located at 

f = 2vre. (80) 
At this fixed point, the theory is a nonrelativistic CFT describing fermions at unitarity. 

2. Scaling dimensions of operators 

The scaling dimensions of the one-fermion operator ip and the two-fermion operator 
(its renormalization is given by Fig. 5) are given by the same formulas as around four 
dimensions; A^, = d/2 and A^^ = 2. Here we concentrate our attention on three-fermion 
operators. The simplest operator is ip^ipi^ip^ which has the orbital angular momentum 
1 = 1. By considering the diagram where one more fermion line is attached to Fig. 5, the 
renormalization of this operator is given by 

/^fp3-^-*W^ + *W^-* 



(27r) d+1 V 2 J V 2 J V 2 

-if It^ ^^ + k-pSc p ^ + k) G ( p -^-k) (81) 



{2n) d+1 V 2 
Zf, N1 A 



(P2 - Pi) In 



47r e- s A 
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FIG. 5: One-loop diagram to renormalize the two-fermion operator ^[ipi- 



Therefore the renormalized operator is (V't'04. VV'iOren = (Z^ ) VtVu^^Tj where 

^t^v*, = 1 + -^s. (82) 

At the fixed point, the anomalous dimension becomes t^^v^t = — <91n Z^^^jds = 
— 3e/2. So the scaling dimension of the operator ■0f0j_V , 0t is 

A^, T ^v^ T = -y + 1 + 7^ T V|V^ T = 4. (83) 

According to the operator-state correspondence, the ground state energy of three fermions 
in a harmonic potential is given by 

[4 + 0(6 2 )]u,. (84) 

For three-fermion operators with I = 0, the calculation is somewhat more involved, be- 
cause there are three lowest operators that can mix with each other: if)^ipi'V 2 ipi, ip^ V^j.- VV>t> 
and ip^ipidtip]. The linear combinations with well-defined anomalous dimensions are 

• V-(V>|'0jVV't) with 7 = -3e/2, 

• ip^il>idtip-\ with 7 = — 3e/2, 

• ^-^j (i9 t + | V 2 ) with 7 = —I. 

The first operator is not a primary operator. Its scaling dimension is trivially equal to 
Ai/'tV'iVi/'t + 1 (corresponding to an excitation of the center of mass motion). The third 
operator annihilates the vacuum and thus does not correspond to any eigenstate of the 
system in a harmonic potential. The second operator is, therefore, the one that corresponds 
to the lowest energy eigenstate of three fermions in a harmonic potential with / = 0. The 
energy of this state is 

u = [5 + 0(e 2 )] u. (85) 



E. 



(0) 



3d „ 3 

h 2 - -e 

2 2 
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We can develop the same analysis for operators having more than three fermion numbers. 
The lowest four-fermion operator with I — is ^^V^-V^. Its anomalous dimension is 
computed to be — 3e, which corresponds to the ground state energy in a harmonic potential 

E ^ = [Q-e + 0(f)]uJ. (86) 

For five fermions, the operator ^•^(VV'fVV'jjV^'-r with the anomalous dimension 
— 7e/2 corresponds to the ground state in a harmonic potential with I = 1. Its energy 
is given by 

= [8-e + 0(f)]uJ. (87) 

We can find two nontrivial operators with / = corresponding to the lowest two energy 
eigenstates of five fermions in a harmonic potential. The operators o^^(Vt(!'|-V^)V 2 )|'f4- 
^tV^i(V^ r V^)V^ t + c^ T ^((ViV^ T )-V^)V^T -#t^i(Vf r Vt/)|)j^ t with 



(a, b, c, d) oc (±19^3 - 5>/35, ^16^/3, -6^35 q= 6>/3, 16>/35) (88) 

have well-defined anomalous dimensions — (51 ± V105) e/16. Therefore, there are five- 
fermion states with / = and energies equal to 



£f = 



lb 



w. (89) 



Finally, the lowest six-fermion operator with / = is ipiipi r (Vipy'Vipi)('Vipy'Vipi). I ts 
anomalous dimension is computed to be — 5e, which corresponds to the ground state energy 
in a harmonic potential 

Ef ] = [l0-2e + O(e 2 )]uj. (90) 

We note that the leading-order results for can be easily understood by recalling that, 
in the limit of d — > 2 from above, fermions at unitarity become noninteracting [10, 12]. So 
the energy eigenvalue of each iV-fermion state is just a sum of single particle energies in a 
harmonic potential at d — 2. Clearly, the ground state energy shows the shell structure at 
d = 2. 



C. Interpolations to d = 3 and discussion 

We determined the exact scaling dimensions of the one-fermion operator ip and the two- 
fermion operator = V'j.V't * n arbitrary spatial dimension d. For the scaling dimensions of 
iV-fermion operators with N > 3, a few lowest-order terms in the expansions over e = d — 2 
and e = 4 — e were computed as summarized in Table I. According to the operator-state 
correspondence, we find that the ground state of three fermions in a harmonic potential has 
the orbital angular momentum / = 1 near d = 2, while I = near d — 4. So there must 
be at least one level crossing between the states with I = and I — 1 as d increases. Using 
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TABLE I: Scaling dimensions of operators corresponding to iV-fermion states with orbital angular 
momentum I in a harmonic potential. Known results for the energy eigenvalues E$ at d = 3 are 
also shown in units of huj [13-15]. 





e = (i — 2 expansion 


e = 4 — d expansion 


known value at d = 3 


N = 2 (1 = 0) 


2 


2 


2 


N = 3 (/ = 0) 


5 + 0(e 2 ) 


4+|e + 0(e 2 ) 


4.66622 


N = 3 (/ = 1) 


4 + 0(e 2 ) 


5-§ e + 0(e 2 ) 


4.27272 


N = 4 (1 = 0) 


6-e±0(e 2 ) 


4 + 8e 2 lnf| + 0(e 3 ) 


5.1 ±0.1 / 5.07 ±0.01 


N = 5 (/ = 0) 


9 11 +f 5I e + 0(e 2 ) 


6+fe + 0(e 2 ) 




iV = 5 (/ = 1) 


8-e + 0(e 2 ) 


7 _n e + 0(e2) 


7.6 ±0.1 


N = 6 (/ = 0) 


10-2e + O(e 2 ) 


6 + 24e 2 lnf| + 0(e 3 ) 


8.7 ±0.1 / 8.67 ±0.03 



the e expansions, the spatial dimension at which this level crossing occurs can be estimated 
to be d ?s 3.4, which means that the three-fermion ground state at d — 3 has I = 1. The 
same level crossing has to occur for the five-fermion case about d ~ 3.64, which implies the 
five-fermion ground state with I — 1 at d — 3. On the other hand, the ground state of four 
or six fermions in a harmonic potential has zero orbital angular momentum near d = 2 and 
d — 4. Thus the level crossing with higher orbital angular momentum states is unlikely and 
we expect I = for the ground state at d — 3. 

In order to make quantitative discussions, we can use the Pade approximants to inter- 
polate the two expansions around d = 2 and d — 4. For each operator, we approximate its 
scaling dimension as a function of d = 2 + e by a ratio of two polynomials; 

L 7 J l±6ie±---±6 y e y v 1 

We demand that the series expansions of (91) around d = 2 and d = 4 match the computed 
results. X ± y is fixed by the number of known terms in the two expansions, while there is 
a freedom in distributing the sum between X and Y. 

The different four Pade approximants for the scaling dimension of each three-fermion 
operator are plotted as functions of d in Fig. 6. We find the behaviors of the Pade approxi- 
mants are quite consistent with the exact results both for I = (left panel) and I = 1 (right 
panel). For three- fermions with / = 0, the interpolated results at d — 3 are 

[3/0] =4.71, [2/1] =4.7, 
[1/2] =4.72, [0/3] =4.72. 

We see that all Pade approximants give very close results in a small interval ps 4.71 ± 
0.01. The harmonic oscillator frequency uo was set to 1 again. This is close to the exact 
result 4.66622 at d = 3 [13], while the numbers obtained by the Pade interpolations are 
slight overestimates of the exact value. 
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FIG. 6: Energy of three fermions in a harmonic potential with I = (left panel) and I = 1 (right 
panel) as functions of spatial dimension d. The four solid curves, although hard to distinguish, 
represent different Pade approximants interpolating the expansions around d = 2 and d = 4. The 
dashed (dotted) lines are extrapolations from the e = 4 — d (e = d — 2) expansions. The symbols 
(x) indicate the exact values for each d [9]. 



Similarly, for three fermions with 1 = 1, the interpolated results are given by 

[3/0] = 4.29, [2/1] = 4.3, 
[1/2] =4.32, [0/3] =4.29, 

which span a small interval E^ w 4.30 ±0.02. A gain the result is very close to, but slightly 
larger than the exact value 4.27272 at d = 3 [13]. One may expect that these agreements 
will be further improved once additional terms in the expansions around d = 2 and d = 4 
are included. 

For four fermions with 1 = 0, the Pade interpolations to d = 3 give 
[4/0] = 5.55, [3/1] = 4.94, [2/2] = 4.94, 

(94) 

[1/3] =4.90, [0/4] = 6.17. 

Some of these estimates are not too far from the numerical result 5.1 ±0.1 [14] and 5.07 ± 
0.01 [15] at d = 3, but the results span a larger interval pa 5.53 ± 0.64. This may 
be because the coefficient of the next-to-leading-order correction around d = 4 is sizable 
compared to the leading term. If one excluded the two extremely asymmetric cases [4/0] 
and [0/4] where all terms in the Pade approximant come to the numerator or denominator, 



(0) 
4 

For five fermions with I = 0, the interpolated results at d = 3 are given by 



one would have a rather small interval about E\ pa 4.92 ± 0.02 



[3/0] = 7.71, [2/1] = 7.64, 
[1/2] = 7.66, [0/3] = 7.82, 



(95) 
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which are in an interval rs 7.73 ± 0.09. On the other hand, for five fermions with / = 1, 
the interpolated results are 

[3/0] = 7.10, [2/1] = 7.16, 
[1/2] = 7.19, [0/3] = 7.09, 

which are in an interval « 7.14 ±0.05. We find E^ < Ef ] at d = 3, which means I = 1 
for the five-fermion ground state in a harmonic potential. However, its energy eigenvalue is 
a substantial underestimate of the numerical result 7.6 ± 0.1 at d = 3 [14]. 
For six fermions with 1 = 0, the Pade interpolations to d = 3 give 

[4/0] = 10.1, [3/1] = 7.92, [2/2] = 7.92, 
[1/3] = 7.80, [0/4] = 16.4. 

The results now span a considerably larger interval E^ m 12.1 ± 4.3, probably because 
of the huge next-to-leading-order coefficient around d = 4. The large error signals the 
worse convergence of the series expansions as the number of fermions increases. If the two 
extremely asymmetric cases [4/0] and [0/4] were excluded, one would have Ef ] m 7.86±0.06. 
For comparison, the numerical result is 8.7 ±0.1 [14] and 8.67 ± 0.03 [15] at d — 3. 

As the number of fermions goes to infinity, one can expect the series expansions over 
e = d — 2 and e = 4 — e for the ground state energy break down. Indeed, the energy of iV 
fermions at unitarity in a harmonic potential scales with different powers of N in different 
spatial dimensions as E N ~ N^ d+V) / d for sufficiently large N. Therefore it is not surprising 
that our extrapolations to d = 3 do not work well for five and six fermions. It is possible 
that the situation is improved once we know the next terms in the expansions around d — 2 
and d = 4. 

Here we comment on the convergence of the e and e expansions. Since the exact integral 
equation to determine the energy eigenvalues of three fermions in a harmonic potential is 
known [9], one can estimate the radii of convergence of the expansions around d = 2 and 
d = 4 by studying their asymptotic behaviors. It turns out that the expansions for the 
three-fermion state with / = are convergent when |e| < 1.0 or |e| < 0.48, while those 
with I = 1 are convergent when |e| < 1.0 or |e| < 1.4. On this basis, we speculate that the 
expansions over e = d — 2 and e = 4 — d have nonzero radii of convergence for systems with 
a finite number of particles. The full details will be reported elsewhere [16] . 



(96) 



IV. EXAMPLE 2: ANYONS 

Anyons in two spatial dimensions present another example of a nonrelativistic CFT. In 
this section, we compute the scaling dimensions of some operators near the bosonic limit 
and the fermionic limit, where perturbative expansions in terms of statistical parameter 9 
are available. Our analytical results, as we will see, are consistent with results obtained 
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FIG. 7: One- loop diagrams to renormalize the contact interaction coupling v. 



by the conventional Rayleigh-Schrddinger perturbation theory in a harmonic potential or 
numerical simulations. 

The field-theoretical representation of anyons is provided by the following Lagrangian 
density where a nonrelativistic field <p is minimally coupled to a Chern-Simons gauge field 

= (a ,a): 

1 1 1 2 

C = Ja dta xa ~ ^ a ° V x a - 77? ( v • a ) 

A9 26 2£ (gg) 

1 2 v 2 

+ i(p* (d t + ia ) <f - - ia) <f\ - - (<p*<p) . 

if is either a bosonic or fermionic field. We denote the propagator of tp by G(p). In the 
Coulomb gauge £ = 0, the only nonvanishing components of the propagator are 

D i0 {p) = -DM = -2x6^. (99) 

We define the three-point vertex T = —1, Ti(p,p') = (p^ + p' t ) /2m and the four-point vertex 
Tij = —5ij/m. The contact interaction coupling v has to be fine-tuned so that the system 
is scale invariant. We start with the case where ip is bosonic. 



A. Near the bosonic limit 



1. Fixed points 



There are two one-loop diagrams which are logarithmically divergent and renormalize the 
coupling v (Fig. 7). Integrating out modes in the momentum shell e~ s A < k < A, the first 
diagram is evaluated as 

r dk AO 2 A 
-4 / G(k)T D 0i (k)^vD j0 (k)T = -i In — — , (100) 

while the second one as 
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FIG. 8: One-loop diagrams to renormalize the two-anyon operators. 



Therefore the renormalization group equation for v is 

dv 4 9 2 v 2 



(102) 



ds n An 
We find two fixed points located at [17] 

v = ±4|0|. (103) 

At these fixed points, the theory is a nonrelativistic CFT. The repulsive (upper sign) or 
attractive (lower sign) contact interaction corresponds to a different boundary condition 
imposed on the s-wave two-body wave function at origin ~ r ± ' 6 '' //7r [18, 19]. 



2. Scaling dimensions of operators 

Since the one-body operator ip is not renormalized, its scaling dimension is A^ = 1, 
independent of 9. Using the operator-state correspondence, we find that there is a one- 
anyon state in a harmonic potential with energy lu, which is obvious. 

The two-body operator ip 2 is renormalized by the diagrams depicted in Fig. 8, which are 
potentially logarithmically divergent and contribute to the anomalous dimension to order 
9. Since the first diagram turns out to be finite, we can concentrate on the second diagram 
which is given by 

i lJ^ k ^ = -^A (104) 

Therefore the renormalized operator differs from the bare operator by a rescaling factor: 
(y^rcn = Z~ 2 V 2 , where 

V = 1 - ^-s. (105) 

47T 

At the fixed point, 7^2 = —dhiZ^/ds = ±\9\/n is the anomalous dimension of the com- 
posite operator ip 2 . We thus find the scaling dimension 

A V 2 =2A ¥ , + 7 ¥ , 2 = 2±M. (106) 
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According to the operator-state correspondence, Lp 2 corresponds to a two-anyon state in a 
harmonic potential with energy equal to 

uo. (107) 

It is straightforward to generalize our analysis to the lowest iV-anyon operator tp N . Its 
scaling dimension is given by 

A ,„ = iV± ^>M, (108) 

Z 71 

and therefore, the corresponding iV-anyon state in a harmonic potential has the energy 
En = A^jv uo. This result coincides with the exact energy eigenvalues of N anyons in a 
harmonic potential [20-22]. 



Eo — 



2±M 



7T 



B. Near the fermionic limit 

1. Scaling dimensions of operators 

If ip is a fermionic field, the contact interaction term (y?*y?) 2 vanishes and the left diagram 
in Fig. 7 turns out to be finite. Therefore, the system is automatically scale invariant. We 
denote the statistical parameter in this case by 9' = 9 — n. 

Let us first look at the lowest two-body operators ipVif. Inserting this operator into the 
left diagram of Fig. 8, its renormalization is given by 



i 



/dk 
(pa - k) i G( Pl + fc)r„(pi,pi + k)D^(k)T 1/ (p 2 ,p 2 - k)G(p 2 - k) 

/dk 
— ( Pl + k\G{ Pl + k)T,( Pl , Pl + k)D, u (k)T u (p 2 ,p 2 - k)G(p 2 - k) (109) 

9' A 
= i-e iJ -(p 2 -p 1 ) j lii— . 

We thus find the linear combinations y?V x ip^itpV y tp with well-defined anomalous dimensions 
7 = ±9'/ir. Therefore, the scaling dimensions of such operators are 

9* 

&<pv x <p^i<pv v <p — 3 ± — ■ (110) 

The operator-state correspondence tells us that the two-anyon states in a harmonic potential 
have the energies 

an 

oo. (Ill) 



Eo 



This result coincides with the exact energy eigenvalues of two anyons in a harmonic potential. 
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Similarly, we can find the lowest three-body operator ipV x ipV y ip has a vanishing anoma- 
lous dimension to order 9'. Therefore, the ground state energy of three anyons in a harmonic 
potential near the fermionic limit is given by 

E 3 = [5 + 0(9' 2 )] to. (112) 

The same result has been derived using the conventional Rayleigh-Schrodinger perturbation 
theory up to order 9' 2 [20]. 

We now turn to the four-anyon case. There are four lowest operators that can mix with 
each other: (pV x (pV y (pV xx (p, ipV x (pV y ipV xy (p , (pV x (pV y (pV yy (p, and (pV x (pV y (pd t ip. The 
linear combinations with well-defined anomalous dimensions are 

• ipV x tpV y tpd t tp with 7 = 0, 

• ip'V x ip'V y ip'V 2 (p with 7 = 0, 

• ipV x ipV y ip (V x - iV y f if with 7 = §£, 

• ipV x (pV y (p (V x + iVyf tp with 7 = ~|f ■ 

The combination of the first two operators y?V x (pV y (p (id t + |V 2 ) ip annihilates the vacuum 
and thus does not correspond to any eigenstate of the system in a harmonic potential. The 
other three operators, therefore, correspond to the energy eigenstates of four anyons in a 
harmonic potential. According to the operator-state correspondence, the energies of these 
three states are given by 

E 4 = [8 + 0(9' 2 )] to and E 4 = 

To our knowledge, these analytical results have not been derived so far. Our numbers are 
consistent with slopes observed in the numerical simulation [23]. 

V. CONCLUSION 

In this paper we study Schrodinger algebra and its representation in terms of operators. 
We show that irreducible representations are built upon primary operators. We also point 
out a correspondence between primary operators and eigenstates in a harmonic potential. 
We illustrate this connection by computing the energy eigenvalues of up to six fermions at 
unitarity in a harmonic potential using expansions over 4 — d and d — 2, as well as the energy 
eigenvalues of up to four anyons in a harmonic potential using expansions over 9 and 9 — ir. 
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±l~ + 0(9' 2 ) 

Z IT 
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TABLE II: Part of the Schrodinger algebra. Given are the values of [A, B\. 



\ I > 

A \ 

\ 


Pi 




D 


C 


H 


Pi 





-iSijN 


-iPi 


-iKi 





Ki 


iSijN 





iKi 





iPi 


D 


iPj 


-iKj 





-2iC 


2iH 


C 


iKj 





2iC 





iD 


H 





-iPj 


-2iH 


-iD 






APPENDIX A: COMMUTATORS IN THE SCHRODINGER ALGEBRA 



The Schrodinger algebra is formed from the operators N, D, M^, Ki, P i: C, H. The 
commutators of N and Mij with other operators are 

[N, D] = [N, M tJ ] = [N, K t ] = [N, P,] = [N, C] = [N, H] = 0, (Al) 



[Ma, M kl ] = i(5 lk M jk - S jk M u + 6 u M kj - 5 3l M ki ), (A2) 

[M^, K k ] = i{5 ik Kj - 8 jk Ki), [M^, P k ] = i(5 tk Pj - 5 jk Pi), (A3) 

[M^, C] = [M^, D] = [M^, H] = 0. (A4) 

The rest of the algebra is summarized in Table II. 



APPENDIX B: THE STRESS TENSOR 



We will find the stress tensor LT^ which appears in the momentum conservation equation: 

d t ji + djllij = 0. (Bl) 
From the evolution equation for the field operators ip and ift it follows that 

d t (j(x)) = i(^t V 2 V^ + V 2 V^^ - VVW - V^ t V 2 ^)(x) 

- JdyV x V(x-y):n(x)n(y):, (B2) 

where : • • • : denotes normal ordering [i.e., :n(x)n(y) := ^(x^^y^piy^aix)}. It is not 
obvious that the right hand side of Eq. (B2) can be written as the derivative of a stress 
tensor. To do that, let us introduce the Laplace transform ct(Q) of the function rV(r): 

oo 

AnrV(r) = JdQa(Q)e- Qr . (B3) 
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In other words, we write the potential V(r) as a superposition of Yukawa potentials, 

oo 

dQa{Q) 6 —. (B4) 



o 



We also introduce, for each value of Q, an auxiliary field <tq(x), 

/ e -Q\x-y\ 
dy An\x-y\ n{y) - (B5) 

It satisfies the equation 

(-W 2 + Q 2 )a Q (x) = n{x). (B6) 
The stress tensor can now be introduced: 

oo 

+ J dQa(Q) :^-d t a Q d,a Q + 5 f [(Va Q ) 2 + Q 2 <Jq] } : . (B7) 
o 

By using Eq. (B6) it is straightforward to verify that Eq. (Bl) is satisfied. 
Notice that ILj is not unique. For example, one can replace 

n - -> IUj + (didj - %V 2 )$ (B8) 

with any $ without destroying the momentum conservation. 

Let us now show that if V(r) is a short-range potential with infinite scattering length, 
then 

dxU ii (x) = 2H. (B9) 
By using Eq. (B6) and the following property of the Yukawa potential, 

p -Q\x-y\ -Q\x-z\ P -Q\y-z\ 

dx ^ r = o 7\ ( B1 °) 

4n\x — y\ 4n\x — z\ 8nQ 
(which can be shown, e.g., by using the Fourier transforms), we find 

dxU ll (x) = 2T + V+ dQdxdya(Q)Q e ^ :n(x)n(y):, (Bll) 

where T is the kinetic energy and V is the potential energy. Since V(r) is a short-range 
potential giving an infinite scattering length, the low-energy physics does not change when 
one rescales the potential as 

V(r) -> AV(Ar). (B12) 

In particular the Hamiltonian is unchanged under the transformation (B12). Setting A = 
1 + e, e < 1 and expanding H to the linear order in e, we find 



/ 



— Q\x— y\ 

dQ dx dya{Q)Q < —- : n(x)n(y) : = V. (B13) 

07T 
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Therefore, we obtain Eq. (B9). 

One can use this relationship to prove that in the normal phase (above the critical tem- 
perature), the bulk viscosity of a Fermi gas at unitarity is identically zero. Indeed, the bulk 
viscosity is given by the Kubo's formula: 



oo 



C = lim -L fdtfdx e^aiL,^, x), %(0, 0)]), (B14) 



o 



but the integral over x can be taken according to Eq. (B9). Moreover, as ([H, O]) — in 
thermal equilibrium for any operator O, the bulk viscosity is zero. This result was derived 
previously using a different approach [24]. 
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